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and will be instrumental in their practical applications. Furthermore, this study aims 
to explore the relationship between the level of indeterminacy and the resulting ran- 


dom variates. By investigating how varying degrees of indeterminacy impact random 
variates, we gain valuable insights into the dynamics of these distributions under dif- 
ferent uncertainty conditions. Preliminary results suggest that random variates exhibit 
a trend of decreasing as indeterminacy levels increase, shedding light on the intriguing 
interplay between indeterminacy and random variate generation. 
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Introduction 

In recent years, due to the complexity of the systems, it may not be possible to evalu- 
ate the work of the real system. The complete study of the complex systems may 
increase the time and cost. When it is not possible to study the system directly, the 
decision-makers from various fields of science depend on simulation. According to 
Jdid et al. [12] “The simulation depends on the application of the study on systems 
similar to the real systems and then projecting these results if they are appropri- 
ate on the real system. The simulation is based on generating a series of random 
numbers that are subject to a uniform probability distribution” Jdid et al. [11] sug- 
gested using probability distribution to generate random numbers for the required 
system. Random numbers are generated through a well-planned generator using the 
probability distribution. The random number generators are run using the algorithm 
designed keeping the mind the target to study the system. A physical or computer 
method that generates random numbers is called the random number generator. 
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A random number generator generates random numbers using the statistical dis- 
tribution and has no specific pattern. Due to high-speed computers, mathematical 
algorithms are designed and implemented to generate random numbers, see [24]. 
Accept-reject simulation method has been widely used to generate random numbers 
using a suitable underlying distribution. In this simulation method, two uniform 
numbers are generated and applied to calculate the probability using the underlying 
distribution. The ratio of probabilities obtained using the first uniform number and 
the mode are compared with the second uniform number in generating the random 
numbers. Devroye [9] suggested the approach for generating random numbers using 
log-concave probability densities. Hurtado and Barbat [10] used the Monte Carlo 
simulation to generate random numbers. Ridout [19] engaged in the development of 
a random number generation technique employing Laplace transformation. Martino 
and Miguez [16] introduced accept/reject sampling to generate random numbers. 
Stein and Keblis [23] introduced the method of simulation for triangular distribu- 
tion. Wang et al. [25] introduced the accept-reject technique for generating random 
variates. Martino et al. [15] worked on accept-reject simulation and presented an 
algorithm to generate random numbers. Mohazzabi and Connolly [17] proposed the 
simulation criteria for normal distribution. Luengo et al. [14] demonstrated param- 
eter estimation through the utilization of Monte Carlo simulation methods. Schinazi 
[20] proposed the method to generate discrete random numbers. Alvarez et al. [2] 
improved the statistical qualities of a random number generator. 

Recently, neutrosophic statistics have gotten attention due to many applications 
for the data recorded from the complex process or when the data has neutrosophic 
numbers, [22]. Many studies including [3, 7, 8] discussed the efficiency of neutro- 
sophic statistics over classical statistics. Alhabib et al. [1] introduced the funda- 
mentals of neutrosophic methods for generating uniform random numbers. Aslam 
[4] introduced the Weibull distribution under indeterminacy as the generalization 
of classical Weibull distribution. Khan et al. [13] explored the gamma distribution 
within the framework of neutrosophic statistics and applied it to analyze complex 
data. Sherwani et al. [21] presented the concept of work distributions under neu- 
trosophic statistics. By following Aslam [4] and Nayana et al. [18] proposed DUS- 
Weibull distribution under indeterminacy. Aslam [6] presented the truncated 
variable algorithm for generating random variates from the neutrosophic DUS- 
Weibull distribution. Second, [5] introduced novel approaches involving sine—cosine 
and convolution techniques to generate random numbers within the neutrosophic 
framework. 

After a comprehensive review of the existing literature and to the best of our 
understanding, no prior research has been conducted on the utilization of the 
accept-reject simulation method within the framework of neutrosophic statistics. 
This paper aims to address this gap by introducing the accept-reject simulation 
technique tailored for neutrosophic uniform distribution (NUD) and neutrosophic 
Weibull distribution (NWD). The following sections will detail the algorithms 
employed for implementing the accept-reject simulation method and will provide 
comprehensive tables of random variates across a spectrum of uncertainty levels. 
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Preliminaries 
In this section, we will introduce the neutrosophic uniform distribution and neutrosophic 
Weibull distribution briefly. 


Neutrosophic uniform distribution (NUD) 

Consider the neutrosophic uniform random variable denoted as xyz, which exhibits a 
degree of uncertainty J,,,, within the range |I,,,J,,,, |. This variable follows a probability 
distribution function (pdf) characterized by indeterminacy, with parameters ay and by 
representing the uniform distribution, as detailed in the work of Aslam [4]. 


1 1 
JS (nu) = (= = =) So (= _ = Foe Len € eaeey aera is, <%nu S bn 
(1) 
Note that the first part denotes the pdf of the uniform distribution under classical sta- 


tistics and the second part denotes the indeterminate part. The corresponding cumulative 


distribution function (cdf), say F (xyz) under indeterminacy can be given as 


xNu — 4u 
| Lic dane water = ine =o 
| Gea a eT 


(2) 


The probability density function, denoted as f («Ns,,), for the case when ay, = ay = ans 


Fnu) = ( 


XNL — AL ) 
(by — az) 


and by = by = bys, and with the degree of indeterminacy Js, can be expressed as follows: 


1 
(bus — ans) 


f (ansy,) = ( Ja + Ins); InseUzs,Jus],4n < xnu < bn (3) 

The streamlined expression for the cumulative distribution function, denoted as F (xNs,,); 
when ay = ay = ans and by = by = bys under conditions of indeterminacy, can be pre- 
sented as follows: 


XNS — 4NS 

Fans) = ( ) (1 + Ins); Ins€Uzs,lus],4n < xnu < bn (4) 
(bns — ans) 

Note the pdf and cdf under indeterminacy reduces to pdf and cdf under classical statistics 


when no indeterminacy is found that is I,,,, = Ins = 0. 


Neutrosophic weibull distribution 

Aslam [4] introduced the concept of the neutrosophic Weibull distribution (NWD) for the 
first time. Let x,y represent a neutrosophic random variable following the Weibull distri- 
bution under indeterminacy, characterized by the scale parameter a, shape parameter f, 
and the measure of indeterminacy Iw within the range [Jz w,, Iw]. The probability density 
function (pdf) for this indeterminate Weibull distribution, as derived in [4], is expressed as 


follows: 


pond {(2) (2) « (EB) artnet tn 


(5) 
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The modified version of the probability density function, denoted as f (xNsy,), for 
the scenario where x, = xy = xs and the measure of indeterminacy is Iys, can be for- 
mulated as follows: 


past xs B 
f (esp) = 1(6) ey ioe) ha + Ins); Insel; lus] (6) 


The cdf corresponding to the Weibull distribution is represented as F (“NSw)s and it 
can be expressed as follows: 


*S 


B 
FtKigge =i = fe) ql ho} + Ins; Ins€Uzs,Ius] (7) 


The mode of the Weibull distribution is represented as f (*yw’) is given by 


1 


a —1)\28 
f(%nw) = «( =) “OP ieeiietid (8) 


The NWD reduces to the Weibull distribution under classical statistics when 
Ins = Inw = 0. 


Neutrosophic accept-reject simulation method 

We present the neutrosophic accept-reject simulation method in the section. In this 
section, we will present the routine to generate random variate from the NUD and 
NWD. Suppose that uyy ~ Uni ((0, 0], [1,1]) and uyn2 ~ Unz2((0, 0], [1,1]) be two ran- 
dom numbers follow the NUD. Suppose that xns,, and xysy, are neutrosophic random 
variables following the NWD and NUD, respectively. Let xyw and xy be the mode 
of xyw and xy, respectively. The random variate xyz will be obtained from the fol- 


lowing expression. 


uN 
(1 + Ins) 


*NSy = 4ns 4 ( ) rs ans); UN€[U Yu], Ins€Uzs, lus] (9) 


The random variate xyw will be obtained from the following expression. 


1 
1— (uy -[ B 
XNSw =o | in( EN “s)) sIns€ Us, lus], un€[4z, Yu] (10) 
1+ Ins 

Based on the above-mentioned information, the following accept-reject routine under 
indeterminacy will be implemented to generate n random variate xys,, from NUD. 

Step-1: pre-specified Ins, ans and bys. 

Step-2: Generate two NUD random numbers, say uni ~ Uni((0,0],[1,1]) and 

un2 ~ Un2((0, 0], [1, 1). 

Step-3: Find xns,, using xys,, = 4ns + (as) (bns — ans) 


Step-4: Compute f («ys,,) using f (xxs,) = ome) (1 + Ins). 


(bns—ans) 
Step-5: Compute f (Ensy) at any value between ays and bys. 
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Step-6: If un2 <f (xns,,)/f (*ns,), then record xy and go to step-7, otherwise, 
repeat the previous steps. 

Step-7: Return «ys, and repeat the process to generate n variate of xys,,. 

The algorithm to generate random variate from the NUD is depicted in Fig. 1. 

The following accept-reject routine will be implemented to generate 1 random vari- 
ate xnsy, from NWD. 

Step-1: pre-specified Ins, ans and bys. 

Step-2: Generate two NUD random numbers, say uy ~ Uni((0,0],[1,1]) and 
ux2 ~ Un2((0,0}, (1, 1). 

Step-3: Find xnysy, using «sy = a [=In( =u) ) . 


ek 
Step-4: Compute f (xysy) using f («nsy) = {(2) (yee (*) ha + Ins). 


a 


Step-5: Compute mode f (%Nsy,) and compute f (nsw) = (4) (ze) \a oe 


Step-6: If un2 <f (xNnsw)/f (%Nsy), then record xyw and go to step-7, otherwise, 
repeat the previous steps. 

Step-7: Return xysy, and repeat the process to generate n variate of xnysyy. 

The algorithm to generate random variate from the NWD is depicted in Fig. 2. 


Note that the proposed neutrosophic accept-reject algorithms are an extension of 


accept-reject algorithms under classical statistics. The proposed simulation methods 
under neutrosophic statistics are reduced to the simulation methods under classical 
statistics when no indeterminacy is found in the data. 


Generate 
Uyi~Uni ([0,0], [1,1]) 
Lyo~ ({0,0], [1,1]) 


Put the values of 
Tvs» Ans and bys 


Start Loop 


Fig. 1 Accept-reject routine for NUD 
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Generate 
un1i~Un1([0,0], [1,1]) 
uUn2~Un2([0,0], [1,1] 


Put the values of 
Ins, Ans and bys 


Start Loop 
Cee) T 


1 
see, <a) in (2a Mtn = Ins 9 
Nsw 1+ Iys 


f (xnsw) 


“Gea "ore 


f(Bnsw) = {(2) (2 (3 (1 + Ins) 


Record Value of 
XNSw 


Fig. 2 Accept-reject routine for NWD 


Simulation studies 

An extensive simulation study has been carried out to investigate the behavior of ran- 
dom variates generated by the developed algorithms across a range of parameters for 
NUD and NWD. In this section, we will present and analyze the simulation tables 
obtained through the accept-reject simulation method and the repetitive simulation 
method. A wide simulation study is conducted to see the behavior of random variate 
obtained from the developed algorithms at various values of parameters of NUD and 
NWD. In this section, the simulation tables obtained from the accept-reject simulation 
method and repetitive simulation method will be given and discussed. 


Accept-reject simulation method 

The random variate from the accept-reject method for the NUD is obtained for vari- 
ous values of ays and bys is reported in Tables 1, 2. To see the effect of intermediacy on 
these random variate, various values of [ys are considered in the generation of random 
variates from NUD and NWD. The random variate obtained from NUD using the values 
of un1, UN2, 4ns=10 and bys=20 are reported in Table 1. The random variate obtained 
from NUD using the values of uy, Un2, ans=15 and bys=20 are reported in Table 2. 
From Tables 1, 2, it can be noted that for the fixed values of ays=15, bys=20 and Ins, 
there is an increase in random variate. The random variate decreases as the values of Ins 
increases from 0 to 1. Note that the values of random variate J)j;=0 presents the values 
of random variate for classical statistics. From Tables 1, 2, it can be noted that inde- 
terminacy plays a significant role in decreasing the random variate. It’s worth empha- 
sizing that computer-based uniform random variate generation is commonly employed 
to simulate random numbers. However, this simulation approach often overlooks the 
influence of the measure of indeterminacy when generating random numbers in the 
computer. Our analysis underscores that, in the presence of uncertainty, the resulting 
random numbers differ from those generated in a deterministic setting. Therefore, when 
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Table 1 Random numbers from NUD when ays=10 and bys=20 


uni =UN2~—s Ins 


0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
0.11 0.28 1.07 0.97 0.89 0.82 0.76 0.71 0.67 0.63 0.59 0.56 0.53 
030 0.50 2.98 271 248 2.29 2.13 1.98 86 1.75 1.65 1.57 1.49 
035 0.59 3.48 3.16 2.90 2.67 248 2.32 2.17 2.05 1.93 1.83 1.74 
0.78 0.56 7.76 7.06 6.47 5.97 5.54 5.17 4.85 457 431 4.09 3.88 
0.14 0.28 141 1.28 1.18 09 01 0.94 0.88 0.83 0.78 0.74 0.71 


085 0.33 8.50 773 7.08 6.54 6.07 5.67 5.3 5.00 472 447 4.25 
0.67 0.52 6.65 6.05 5.54 5.12 4.75 443 4.16 3.91 3.70 3.50 3.33 
0.15 0.20 1.55 141 1.29 19 10 1.03 0.97 0.91 0.86 0.8 0.77 
0.18 048 1.78 1.62 1.49 37 .2/ 1.19 J 1.05 0.99 0.94 0.89 
087 061 8.69 7.90 7.24 6.69 6.21 5.80 5.43 5.11 483 458 435 
0.76 0.02 7.62 6.92 6.35 5.86 5.44 5.08 476 448 4.23 40 3.81 
0.61 0.75 6.09 5.54 5.08 4.69 435 4.06 3.8 3.58 3.38 3.2 3.05 
0.67 0.28 6.73 6.12 5.6 5.18 481 449 42 3.96 3.74 3.54 3.37 
0.06 0.88 0.55 0.50 0.46 0.43 0.40 0.37 0.35 0.33 0.31 0.29 0.28 
044 0.92 438 3.98 3.65 3.37 3.13 2.92 2.74 2.58 243 23 2.19 


unit =-UN2~—s INS 


0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 


0.11 0.28 5.53 5.49 5.44 541 5.38 5.36 5.33 5.31 5.30 5.28 5.27 
030 0.50 6.49 6.35 6.24 6.14 6.06 5.99 5.93 5.88 5.83 5.78 5.74 
035 0.59 6.74 6.58 6.45 6.34 6.24 6.16 6.09 6.02 5.97 5.92 5.87 
0.78 0.56 8.88 8.53 8.23 7.99 777 759 743 7.28 7.16 7.04 6.94 
0.14 0.28 5.71 5.64 5.59 5.54 5.50 547 5.44 541 5.39 5.37 5.35 
085 0.33 9.25 8.86 8.54 8.27 8.04 783 7.66 7.50 7.36 7.24 7.12 
0.67 0.52 8.33 8.02 7.77 7.56 7.38 7.22 7.08 6.96 6.85 6.75 6.66 
0.15 0.20 5.77 5.70 5.64 5.59 5.55 5.52 5.48 5.45 5.43 541 5.39 
0.18 048 5.89 5.81 5.74 5.69 5.64 5.59 5.56 5.52 5.50 5.47 5.45 
087 061 9.35 8.95 8.62 8.34 8.10 7.90 772 7.56 7A\ 7.29 7.17 
0.76 0.02 8.81 8.46 8.17 7.93 772 7.54 7.38 7.24 7.12 7.00 6.90 
0.61 0.75 8.05 7.77 7.54 7.34 7.18 7.03 6.90 6.79 6.69 6.60 6.52 
0.67 0.28 8.37 8.06 78 7.59 7A) 7.24 7.10 6.98 6.87 6.77 6.68 
0.06 0.88 5.28 5.25 5.23 5.21 5.20 5.18 5.17 5.16 5.15 5.15 5.14 
044 0.92 7.19 6.99 6.83 6.68 6.56 6.46 6.37 6.29 6.22 6.15 6.10 


faced with uncertainty, depending on uniformly distributed random variables generated 
using traditional simulation methods could result in erroneous conclusions. 

The random variate from NWD for various values of Js and parameters are shown 
in Tables 3, 4, 5, 6. Table 3 is presented for a = 2 and 6 = 0.50. The random variate for 
a =2 and $ = 1 (exponential distribution) is shown in Table 4. The random variate 
for @ = 2 and 8B = 1.5 is shown in Table 5. The random variate for a = 2 and 6B = 1.5 
is shown in Table 5. The random variate for a = 2 and B = 2 is shown in Table 5. 
From Tables 3, 4, 5, 6, there is a decreasing trend in random variate as the values of 
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Table 3 Random numbers from NWD when a@ = 2 and B = 0.50 


uN1 


UN2 


Ins 


0.1 


0.2 


0.3 


0.4 


0.5 


0.7 


0.8 


0.1268 
0.5698 
0.7631 
0.1133 
0.6222 
0.5988 
0.9762 
0.0749 
0.6036 
0.7040 


0.4469 
0.6751 
0.8754 
0.6995 
0.2475 
0.7014 
0.5398 
0.6296 
0.7176 
0.2784 


0.0300 
0650 
0.9224 
0.0236 
3910 
2361 
0.9224 
0.0100 
2059 
2.0876 


0.0249 
0.8294 
0.7230 
0.0197 
1.0685 
0.9556 
0.7230 
0.0083 
0.9775 
1.5612 


0.0211 
0.6653 
0.5827 
0.0166 
0.8484 
0.7623 
0.5827 
0.0071 
0.7791 
1.2165 


0.0180 
0.5460 
0.4800 
0.0142 
0.6910 
0.6231 
0.4800 
0.0061 
0.6364 
0.9769 


0.0156 
0.4565 
0.4025 
0.0123 
0.5743 
0.5193 
0.4025 
0.0053 
0.5300 
0.8030 


0.0120 
0.3332 
0.2951 
0.0095 
0.4154 
0.3772 
0.2951 
0.0041 
0.3847 
0.5717 


0.0107 
0.2897 
0.2570 
0.0085 
0.3598 
0.3272 
0.2570 
0.0036 
0.3336 
0.4923 


0.0095 
0.2542 
0.2258 
0.0076 
0.3148 
0.2867 
0.2258 
0.0032 
0.2922 
0.4285 


Table 4 Random numbers from NWD when a = 2and B = 


uN 


UN2 


Ins 


0 


0.1 


0.2 


0.3 


0.4 


0.5 


0.7 


0.8 


0.9 


0.6724 
0.3384 
0.9633 
0.0344 
0.7971 
0.8502 
0.7074 
0.3028 
0.9702 
0.8237 


0.393 
0.4095 
0.8596 
0.1764 
0.8033 
0.750 
0.7036 
0.3204 
0.5936 
0.2993 


0.9224 
0.7352 
0.9224 
0.0635 
0.9224 
0.9224 
0.9224 
0.6440 
0.9224 
0.9224 


0.7230 
0.6625 
0.7230 
0.0581 
0.7230 
0.7230 
0.7230 
0.5816 
0.7230 
0.7230 


0.5827 
0.6030 
0.5827 
0.0536 
0.5827 
0.5827 
0.5827 
0.5304 
0.5827 
0.5827 


0.4800 
0.5533 
0.4800 
0.0497 
0.4800 
0.4800 
0.4800 
0.4874 
0.4800 
0.4800 


0.4025 
0.5113 
0.4025 
0.0463 
0.4025 
0.4025 
0.4025 
0.4510 
0.4025 
0.4025 


0.34 


25 


0.4752 


0.34 
0.04 
0.34 
0.34 
0.34 
0.4 
0.34 
0.34 


25 
34 
25 
25 
25 
96 
25 
25 


0.295 
0.4439 
0.295 
0.0408 
0.295 
0.295 
0.295 
0.3923 
0.295 
0.295 


0.2570 
0.4165 
0.2570 
0.0385 
0.2570 
0.2570 
0.2570 
0.3684 
0.2570 
0.2570 


0.2258 
0.3922 
0.2258 
0.0365 
0.2258 
0.2258 
0.2258 
0.3472 
0.2258 
0.2258 


Table 5 Random numbers from NWD when @ = 2 and B = 1.50 


uN1 


UN2 


Ins 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


0.1946 
0.1403 
0.1072 
0.8911 
0.8834 
0.2193 
0.3930 
0.4797 
0.2676 
0.6454 


0.5075 
0.2019 
0.6891 
0.8086 
0.0830 
0.1400 
0.5098 
0.6271 
0.1190 
0.1619 


0.6717 
0.5301 
0.4380 
0.9224 
2.7644 
0.7339 
1.1606 
1.3794 
0.8533 
1.8417 


0.6303 
0.4982 
04121 
0.7230 
24217 
0.6880 
1.0800 
1.2772 
0.7985 
1.6829 


0.5947 
0.4708 
0.3897 
0.5827 
2.1800 
0.6487 
1.0121 
1.1925 
0.7517 
1.5558 


0.5637 
0.4468 
0.3702 
0.4800 
1.9959 
0.6145 
0.9542 
1.1207 
0.7113 
1.4512 


0.5365 
0.4257 
0.3529 
0.4025 
1.8491 
0.5846 
0.9039 
1.0590 
0.6759 
1.3630 


0.51 


23 


0.4069 


0.33 


75 


0.3425 
1.7282 


0.55 
0.85 


80 
98 


1.0053 
0.6446 
1.2875 


0.4907 
0.3901 
0.3237 
0.2951 
1.6263 
0.5343 
0.8208 
0.9580 
0.6168 
1.2218 


0.4713 
0.3749 
0.3112 
0.2570 
1.5388 
0.5130 
0.7860 
0.9159 
0.5917 
1.1641 


0.4537 
0.3611 
0.2999 
0.2258 
1.4626 
0.4936 
0.7546 
0.8782 
0.5691 
1.1128 


Ins increase. For example, for the exponential distribution, when uy = 0.6724 and 
un2 = 0.3931, the value of random variate is 1.2212 when Jys=0 and the value of ran- 
dom variate is 2.000. 
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Table 6 Random numbers from NWD when @ = 2 and B = 2.0 


UNI un2 Ins 


0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 


0.9215 0.6408 1.2212 0.9224 0.7230 0.5827 04800 04025 0.3425 0.2951 0.2570 0.2258 0.2000 
0.7373 0.0830 2.3125 2.1068 1.9525 1.8303 1.7297 1.6449 1.5719 1.5082 14519 14016 1.3564 
04176 0.5260 14705 1.3820 1.3080 1.2450 1.1904 1.1424 1.0999 1.0619 1.0276 0.9964 0.9679 
0.5954 0.3974 1.9025 1.7655 1.6559 1.5652 14885 14223 1.3644 1.3132 1.2675 1.2263 1.1889 
0.2674 0.0462 1.1155 1.0554 1.0041 0.9597 0.9207 0.8861 0.8552 0.8273 0.8020 0.7788 0.7576 
0.0808 0.4036 0.5806 0.5525 0.5281 0.5067 04877 04707 0.4554 04414 04287 04170 0.4062 
0.5406 0.9562 1.2212 0.9224 0.7230 0.5827 04800 04025 0.3425 0.2951 0.2570 0.2258 0.2000 
04129 0.6503 14594 1.3719 1.2987 1.2363 1.1822 1.1348 1.0926 1.0549 1.0209 0.9899 0.9617 
0.8789 0.6710 1.2212 0.9224 0.7230 0.5827 04800 04025 0.3425 0.2951 0.2570 0.2258 0.2000 
0.7306 0.5675 2.2904 2.0891 1.9376 1.8171 1.7179 1.6341 1.5620 14989 14431 1.3933 1.3485 


Concluding remarks 

This paper has introduced two simulation methods for generating random variates using 
NUD and NWD. The paper presented the essential computational procedures, simula- 
tion techniques, and algorithms necessary for these methods. Random variates derived 
from NUD and NWD were generated under various parameter settings and indetermi- 
nacy levels. The analysis reveals that the measure of indeterminacy significantly influ- 
ences the random variates, with a decrease observed as the level of indeterminacy 
increases. It’s important to highlight that computer-based uniform random variate gen- 
eration is a prevalent method for simulating random numbers. Nonetheless, this simula- 
tion approach frequently disregards the impact of the measure of indeterminacy when 
generating random numbers in the computer. Our analysis underscores that in the pres- 
ence of uncertainty, the resulting random numbers deviate from those generated in a 
deterministic environment. Therefore, depending on uniform random variables gen- 
erated through traditional simulation methods can potentially lead to erroneous deci- 
sions. These proposed simulation methods have versatile applications in fields such as 
education, medicine, computer science, and engineering. Furthermore, the potential for 
extending these methods to handle other distributions is a promising avenue for future 
research. 
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